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Abstract: The genetic diversity and population structure of the Dy¬ 
bowski’s frog {Rana dybowskii ) were investigated by using 11 polymor¬ 
phic microsatellite loci. Total 75 individuals were sampled from six 
populations in Lesser Khingan Mountains and Changbai Mountains, 
China. Results showed that allele number of the 11 microsatellite loci 
was in the range of 2-10 in all populations, with the mean of 5.6. The 
average expected heterozygosity (H E ) was 0.572, indicating a moderate 
polymorphism. The results of genetic differentiation coefficient (Fst) 
showed that population genetic differentiation was significant between 
Changbai and Lesser Khingan Mountains (p<0.001). This result was 
verified further by Nei’s genetic distance (D A ) based on UPGMA phy¬ 
logenetic trees and by AMOVA analysis. In conclusion, the populations 
distributed in Lesser Khingan Mountains and Changbai Mountain are 
proposed to be two distinct management units (MUs) for their protection 
and management. 

Keywords: Dybowski’s frog {Rana dybowskii ); genetic diversity; genetic 
differentiation; management units (MUs); microsatellite 

Introduction 

Dybowski’s frog (Rana dybowskii), is an important economic 
species distributed in Northeast China (Lesser Khingan Moun¬ 
tains and Changbai Mountain), Far East of Russia, Korean pen¬ 
insula, and Japan (Xie et al. 1999). The taxonomy of this species 
was disputable in history. It was identified as a subspecies of 
Chinese brown frog (Rana chensinensis) before 1999 (Boring 
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1945; Orlova et al. 1977; Pope 1931; Stejneger 1925; Wei et al. 
1991a). After then, it was proposed to be an independent species 
according to its difference of distribution range and morphology 
from the R. chensinensis (Xie et al. 1999). This proposal was 
further proved by the analysis of mitochondria DNA (mtDNA) 
control region that addressed the significant genetic difference 
between the two species (Jiang and Zhou 2001; Jiang et al. 2002; 
Yang et al. 2001). 

Meanwhile, it was observed in accumulating studies that the 
Dybowski’s frog in different regions of Northeast China has 
variations in morphology, cytology, and genetics. The female 
frogs in Changbai Mountains have a high reproductive capacity 
and oviduct mass (The oviduct of the frog is described as frog oil 
in traditional Chinese medicine) (Li et al. 2003). The position of 
anus of the tadpole in Liaoning lies closer to its head relative to 
that in Changbai Mountains (Wang et al. 2006). Body size of the 
adult frogs (both male and female) in Lesser Khingan Mountains 
is larger relative to the frogs in Changbai Mountains and Liaon¬ 
ing regions. The similar trend was also observed in the body size 
and body mass of tadpoles (Table 1). A cluster analysis based on 
a set of morphological characteristics further proved the signifi¬ 
cant variation between populations in Changbai Mountains and 
Lesser Khingan Mountains (Ying et al. 2008). Moreover, cen¬ 
tromere indexes of the eighth chromosome in frogs in Liaoning 
differed largely from that of the other regions (Table 1) (Shao et 
al. 1999). The sequence variation of mitochondrial cytochrome b 
gene was 0.3%—1.8% between Liaoning populations and 
Heilongjiang populations (Yang et al. 2001). However, genetic 
distances amongst eight populations in Heilongjiang regions, 
including Yichun, Shangzhi, Wudihe, Shanhe, Suiyang, 
Fangzheng, Chaihe, and Dongfanghong, inferred from RAPD 
was not completely matched with their actual, geographical dis¬ 
tances (Xiao et al. 2001). Meanwhile, no significant genetic dif¬ 
ferentiation was detected in subpopulations in the central 
Changbai Mountains (Li et al. 2009). These results indicated 
significant differentiation among populations in Lesser Khingan 
Mountains, Changbai Mountains and Liaoning region, but not 
within these three regions. However, genetic markers used in 
these studies are not sensitive and reliable to address genetic 
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differentiation in this case. For example, Cyt b gene of mtDNA 
is commonly used in interspecies studies and not sensitive to 
detect subpopulation differentiations (Yang et al. 2001). RAPD 
and RAMP heavily rely on template DNA quality that is not 
applicable to wild populations due to difficulties in sampling and 
sample preservation prior to lab experiments, and now nearly is 
given up in the analysis of genetic diversity and differentiation. 
Therefore, genetic diversity and genetic structure of the Dy- 
bowski’s frogs in different regions of Northeast China requires 
further study by using sensitive and reliable molecular markers. 


Microsatellite is a kind of short tandom repeat DNA in the 
genome with high polymorphism (Sota and Sasabe 2006). PCR 
based analysis is robust and endurable for low quality samples 
(Flolst-Jensen et al. 2003). It has been used in tests of genetic 
diversity and population differentiation in numerous studies and 
its validity for this issue has been proved. Therefore, we carried 
out a systematic study on the genetic diversity and differentiation 
of the Dybowski’s frog in Northeast China by using microsatel¬ 
lites. 


Table 1. Morphological and Chromosomal differences in tadpoles and adult Dybowski’s frog from different geographic origins 


Distribution 

Female 

Average 

Length 

(mm) 

Female 

Average 

Weight 

(g) 

Male 

Average 

Length 

(mm) 

Male 

Average 

Mass (g) 

Length of 
Tadpoles 

(mm) 

Width of 

Tadpoles 

(mm) 

Height at 

Withers of 

Tadpoles 

(mm) 

Tail Length of Hindlimb 

Tadpoles Length of 

(mm) Tadpoles 

(mm) 

CTC 

Heilongjiang 

64.64i2.27 [1] 

31.6i5.59 ™ 

56.7i4.1 121 






M [SI 

Jilin 

>65 [31 

28.3i7.10 [41 

54.2i3.9 [21 

24.70i2.18 m 

34.63±2.97 m 

8.27±0.70 m 

7.39i0.77 [?I 

22.37i2.18 m 0.94i0.05 m 

M 191 

Liaoning 

>63 151 

22.6i2.3 161 

>51 151 

19.4i2.0 161 

20.85i0.93 m 

5.41i0.34 m 

4.48i0.43 m 

12.8i0.75 m 0.63i0.03 171 

St 1101 


[1) (Zhang et al. 2008), PI (Xie et al. 1999), [31 (Wei et al. 1991b), m (Zhao 1982), [5] (Yu et al. 2000), [6) (Tian et al. 2001), m (Wang et al. 2006), [SI (Luo; Li 1985), [91 
(Wei et al. 1991a), [101 (Shao et al. 1999). 

“CTC” means Centromere Type of No. 8 Chromosome “M” means metacentric chromosome; “St” means acrocentric chromosome. 


Materials and methods 

Sample collection and DNA isolation 

A total 75 Dybowski’s frogs were captured in from six popula¬ 
tions distributing in the Songhua River Basin, the Mudanjiang 
River Basin, and the Yalu River Basin in Lesser Kliingan Moun¬ 


tains and Changbai Mountains (Fig. 1). Frogs were killed by 
destroying brain and spinal cord. Muscles were isolated from the 
hind leg of each individual and stored liquid nitrogen immedi¬ 
ately and moved to -80°C refrigerator after arrived at lab. Infor¬ 
mation of samples is shown in Table 2. Genomic DNA was ex¬ 
tracted by standard phenol-chloroform protocol (Elphinstone et 
al. 2003) for all samples, and quantified by using DU-640 nu¬ 
cleic Acid and Protein Analyzer (Beckman USA). 



Fig. I Sampling sites of Dybowski’s frog, (a) Landscape; (b) River system 


Microsatellite amplification 

Eleven microsatellites of genus Rana were selected, including 
RpilOO, RpilOl, Rpil02, Rpil04, Rpil07 (Hoffman et al. 2003), 
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RsyC23, RsyC52, RsyC41, RsyD25, RsyD40, RsyD88 (Julian and 
King 2003) (Table 3). The core sequences of these microsatellite 
loci consist of four nucleotide repeat. The 5’ end of upper prim¬ 
ers of all loci was fluorescently labeled with 5’-FAM, TAMRA, 
and JOE (Shanghai Invitrogen Bio-Technology Co. Ltd). PCR 
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reactions were set up in 10-pL system containing 1 x PCR buffer, 
0.3 mM dNTP, 2.5 mM MgCE, 0.3 pM of each of forward and 
reverse primer, 0.5 units of Taq polymerase (TaKaRa, Dalian, 
China). Cycling was initiated with incubation at 94°C for 5 min, 
followed by 40 cycles of 94°C for 30 s, appropriate annealing 
temperature (Table 3) for 30s, 72°C for 30s, and a final extension 
step of 10 min at 72°C. Then 2 pL of the amplified products were 


mixed with 20 pL of Hi-Di formamide and 1 pL of internal stan¬ 
dard (ILS600, Promega). After denaturation at 95°C for 5 min, 
they were subjected to capillary electrophoresis on ABI 3130 
Genetic Analyzer (Applied Biosystems, USA). Signals were 
collected and processed using software GENEMAPPER 3.2 
(Applied Biosystems, USA). Alleles were defined on the basis of 
fragment size and the length of core sequence. 


Table 2. Summary of information of samples used in this study 


Population sampling site 

Distribution 

River Basin 

Sampling year 

No. of Samples 

Tonghua, Jilin 

Changbai Mountains 

Yalu River Basin 

2006.10 

13 

Jiaohe, Jilin 

Changbai Mountains 

Songhua River and a tributary of Mudan River 

2006.10 

10 

Huadian, Jilin 

Changbai Mountains 

Songhua River Basin 

2008.05 

12 

Mudanjiang, Heilongjiang 

Changbai Mountains 

Mudan River (main tributary of Songhua River) 

2008.05 

10 

Hebei, Heilongjiang 

Xiaoxing’an Mountains 

Songhua River Basin 

2006.10 

15 

Tieli, Heilongjiang 

Xiaoxing’an Mountains 

Songhua River Basin 

2006.10 

15 


Table 3. Primer sequences and annealing temperatures of microsatellite loci used in this study 


Locus 


RpilOO 

RpilOl 

Rpil02 

Rpil04 

Rpil07 

RsyC23 

RsyC41 

RsyC52 

RsyD25 

RsyD40 

RsyD88 


Forward 5'-3' 


CACCCTTAAAAGGACAGAACATT 

C GTTAACGC AC AGCAAAGGAGT A 

GTGTGTGTGTTTATTTACTG 

CCTGATAAAGGGGTTTGTGAAT 

CCGAGGTACCTAGTTGATGTG 

AGGGCATTATTACATTTTGGTC 

GGCAGTCTGGTCCAGTCGTCT 

C CAT AC AACCGTG ATT AC AAAAG 

GACCAGAAAGTTATTTCAAGGG 

TGATTGATTGTTCACTATTGGG 

TCAATCCATCAGTCTGTCTGTC 


Reverse 5'-3' Repeat motif Annealing 

temperature 

_ CQ 

ACCTCTTATTTGTGCCTAACTGAA (GATA) 10 -(GATA) 2 -(GATA) 2 -(GATA)2 63 

GCATGGACAAGGGATGACTTAGAA (GATA) 13 58 

CTTCCATTTTAATTGTGT (GATA),g 50 

GAACCATAAAATGTTGGGATAGAG (AGAC)i 0 -(GATA) 14 50 

CCGAGGTACCGTGATTATGT (GATA) 16 -(GATA) 10 -(GACA) 5 -(GATA) 12 59 

AGGAAATTACAGAGGACTGTGG (TACA) 10 57 

CCACAAAACAGGAATCGGTCATA (TACA)g 62 

ATATACCACCCTTCCAGAGATG (TACA) 17 57 

CCCTGTAACTATGTACCAGGAGG (TAGA) lg 61 

AAGTAGATTATGTGCTGCAAACTG (TAGA) 19 55 

GGATTTTGTAAAGAATGCTCCTC (TAGA) 13 _60 


Data analysis 

Statistical parameters were calculated by using FSTAT 2.9.3 
(Goudet 2001) for 11 microsatellite loci, including number of 
alleles (N A ), observed heterozygosity ( H 0 ), expected heterozy¬ 
gosity ( H e ) and allelic richness (A R ). Fixation index (F IS ) of each 
population was estimated by the combination of FSTAT 2.9.3 
and ARLEQUIN 3.01 (Weir and Cockerham 1984; Excoffier et 
al. 2005). Polymorphism information content ( PIC) per loci was 
estimated based on allelic frequency as described in Xing et al. 
(Xing 2006). The Hardy-Weinberg equilibrium and the het¬ 
erozygote excess or deficiency of each locus in each population 
was tested using GENEPOP 3.4 (Rousset and Raymond 1995; 
Gou and Thompson 1992). Moreover, Linkage disequilibrium 
per locus was analyzed by Fisher's test. 

Genetic differentiation coefficient (Fst) was calculated to 
analyze genetic structure and differentiation of populations by 
using ARLEQUIN 3.01. Significance was tested based on 1,000 
replicates random bootstrap. Interpopulation gene flow (N e m) 
was estimated on the basis of the equation 7V e m=(l — -F ST )/(4.F ST ) 
(Wright 1965). Nei’s genetic distance (Z) A ; Nei et al. 1983) of 


each population was calculated using POPULATIONS 1.2 
(Olivier Langella 2000). Phylogenetic tree was constructed by 
using UPGMA clustering, and its accuracy was estimated by the 
percentage of bootstrapping 1,000 replicates. Genetic differentia¬ 
tion between populations was analyzed by molecular variance 
(AMOVA) to verify obtained clustering results. 

Results and analysis 

Genetic diversity 

A total 62 alleles were observed on 11 microsatellite loci in 75 
samples. Average allele number per loci was 5.6. The allelic 
number on RsyD40 (A a = 10) was the highest among 11 microsa¬ 
tellite loci (Table 4), followed by RsyD88 (N A = 9) and RsyD25 
(N a = 9). RpilOO (N a = 2) had the lowest allelic number. One allele 
unique to “Tieli” population was observed on RsyC23 and 
Rpil04 with frequencies of 0.033 and 0.100, respectively. A 
unique allele with a frequency of 0.133 was observed on RsyD40 
in “Hebei” population, and on locus Rpil04 of “Tonghua” popu¬ 
lation had a unique allele with a frequency of 0.077. 
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* 0.01<P<0.05; **P<0.01 
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The H„ across 11 loci ranged from 0.108 ( RpilOl ) to 0.473 
( Rpil04 ), with an average of 0.291. The H E was between 0.306 
{RpilOl) and 0.827 ( RsyD40 ), averaged 0.572. The averaged A R 
and PIC were 4.475 and 0.561, respectively (Table 4). H„ was the 
highest (H o =0. 336) in “Tonghua” population and lowest (H = 
0.255) in “Mudanjiang” population. However, no significant 
difference for H E was observed among six populations (H E was 
between 0.550 and 0.590). 

Hardy-Weinberg equilibrium and linkage disequilibrium 

As estimated by FSTAT 2.9.3, _F IS , a measure of the deviation 
from random mating within the 6 populations, was 0.504, and 
ranged from -0.301 to 1.000. Among all 11 loci, RpilOO, Rpil07, 
RsyC23 and RsyC52 fit H-W equilibrium (p> 0.05, F ls varied 
from -0.301 to 0.640) in all populations investigated. Locus 
RpilOl, Rpil02 and Rpil04 significantly deviated from H-W 
equilibrium (p<0.01) in “Huadian”, "Hebei”, and “Tieli” popula¬ 
tions, and the rest four loci ( RsyC41 , RsyD25, RsyD40 and 
RsyD88) significantly deviated from the equilibrium (p<0.01, 
average F ls = 0.614, 0.629, 0.667 and 0.775, respectively) in all 
populations (Table 4). 

Totally 55 locus pairs which were randomly combined from 11 
microsatellite loci were tested for linkage disequilibrium. All 
locus pairs but three {Rpil07-Rpil02, RpilOl -RsyC41, 
RsyD25-RsyC41 ) did not show significant linkage disequilibrium 
(p>0.05), accounting for 94.55% of the total pairs. 

Population genetic differentiation 

Genetic differentiation coefficients (Fst) an d gene flows (N e m) 
among six populations were calculated on the basis of the allelic 
frequency of 11 microsatellites (Tables 5 and 6). F ST values 
ranged from 0.040 to 0.178, and A e m values were in the range of 
2.57 and 12.35. Significant genetic differentiation occurred be¬ 
tween the populations in Lesser Khingan Mountains and Chang- 
bai Mountains. Moreover, there was an extremely significant 
genetic differentiation between "Hebei” and "Tieli” populations 
in Lesser Khingan Mountains (p<0.001), whereas genetic differ¬ 
entiation among the populations in Changbai Mountains was not 
significant (p>0.05). 

Genetic distance ( D A ) -based phylogenetic tree showed that 
the populations in Lesser Khingan Mountains, namely "Hebei” 
population and "Tieli” population, fell into one cluster (Table 6 
and Fig. 2), with the bootstrap test value of 92%, and the four 
populations in Changbai Mountains were orderly grouped in a 
cluster, with the bootstrap test value of 85%. This structure was 
accordance with the geographic relationaships among popula¬ 
tions (Fig. 1). 

The six populations were divided into two groups on the basis 
of the above results for AMOVA analysis (Table 7). The results 
showed that genetic differentiation among populations within the 
same group and among individuals within one population were 
both extremely significant (p<0.001). Although there was a small 
genetic variation (/ 7 CT =0.057) between the populations of Lesser 
Khingan Mountains and Changbai mountains, their genetic dif¬ 


ferentiation was still at an extremely significant level (p=0.008). 



Fig. 2 Phylogenetic (UPGMA) tree of six populations based on Nei’s 
genetic distance /), 


Table 5. Pairwise h\i values (below diagonal) and the probability 
(above diagonal) of Dybowski’s frog among six populations 


Population 

TH 

JH 

HD 

MDJ 

HB 

TL 

TH 


** 

*** 

* 

*** 

*** 

JH 

0.082 


NS 

NS 

*** 

*** 

HD 

0.125 

0.041 


NS 

*** 

*** 

MDJ 

0.071 

0.040 

0.062 


*** 

*** 

HB 

0.131 

0.104 

0.178 

0.158 


*** 

TL 

0.129 

0.100 

0.139 

0.132 

0.101 



*** P0.001; ** P<0.01; * P<0.05; NS-not significant 


Table 6. Pairwise estimates of gene flow (Aim, below diagonal) and 
Nei’s D a distance (above diagonal) of the Dybowski’s frog among six 
populations 


Population 

TH 

JH 

HD 

MDJ 

HB 

TL 

TH 


0.144 

0.185 

0.129 

0.184 

0.211 

JH 

5.83 


0.685 

0.097 

0.171 

0.166 

HD 

3.75 

12.01 


0.132 

0.207 

0.183 

MDJ 

6.84 

12.35 

7.79 


0.194 

0.180 

HB 

3.57 

4.55 

2.57 

2.91 


0.103 

TL 

3.64 

4.75 

3.35 

3.55 

4.72 



Table 7. Analysis of molecular variance (AMOVA) 


Source of variation 

d.f. 

Sum of 

squares 

Variance 

components 

Percentage 

of varia¬ 
tion 

Fixation 

Indices 

Among groups 

1 

83.417 

0.201 

5.73 

0.057** 

Among populations 

4 

95.979 

0.263 

7.50 

0.080*** 

within groups 






within populations 

542 

1652.651 

3.049 

86.78 

0.132*** 


***p<0.001; ** 0.001<p<0.01 


Discussions and recommendations 

Genetic diversity 

Genetic diversity is an important genetic parameter for evalua¬ 
tion of evolutionary potential and exploitation potential. There¬ 
fore, it is the theoretical basis of conservation strategies and 
policies. In present study, the average allelic number of 11 mi¬ 
crosatellite loci was 5.6, and the average A R and PIC were 4.475 
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and 0.561, respectively. 

The H e value at species level ranged from 0.306 to 0.827 
across 11 loci, with the mean of 0.572. This figure was signifi¬ 
cantly higher than that in Italian agile frog ( Rana latastei, H E , 
0.163-0.242) (Ficetolal et al. 2007), moor frog ( Rana arvalis, H E , 
0.160-0.552) in Sweden and Denmark (Knopp et al. 2007), relict 
leopard Frog ( Rana onca , H E , 0.11-0.71) (Savage and Jaeger 
2009), wood frog ( Rana sylvatica , H E , 0.14-0.78) (Newman and 
Squire 2001), Columbia spotted frog ( Rana luteiventris, H E , 
0.23-0.70) (Funk et al. 2005), and common frog ( Rana tempo- 
raria, H E , 0.16-0.36) (Matsuba and Merila 2009), and approxi¬ 
mately equal to that in Cascades frog ( Rana cascadae, H E , 
0.25-0.87) (Monsen and Blouin 2003) and moor frog ( Rana ar¬ 
valis, H e , 0.50-0.68) in the Netherlands (Arens et al. 2007), but 
slightly lower than that of Chinese wood frog ( Rana chensinensis, 
H e , 0.504-0.855) (Zhan et al. 2009). This demonstrates that the 
Dybowski’s frog in Northeast China has abundant genetic diver¬ 
sity. 

Population genetic differentiation 

Genetic distance between populations presents genetic variation 
of a species. The phylogenetic tree (Fig. 2) showed the six popu¬ 
lations were clustered to two reliable clusters: one composed of 
"Hcibei” and “Tieli” populations in Lesser Khingan Mountains, 
the other consisted of four populations in Changbai Mountains. 
This demonstrated clear differentiation between populations in 
Lesser Khingan Mountains (bootstrap test value, 92%) and 
Changbai Mountains (bootstrap test value, 85%). AMOVA 
analysis showed that genetic variation caused by the geographic 
difference accounted for 5.73% of total genetic variation. Al¬ 
though the contribution is relatively small, but it caused signifi¬ 
cant genetic differentiation (p=0.008). This suggested that local 
populations had evolved unique genetic characteristics, espe¬ 
cially between Lesser Kingan Mountains and Changbai Moun¬ 
tains. The strength of gene flow was correlation with the geo¬ 
graphic distance among populations. As shown in Table 6, the 
least gene flow (7V e m=2.57) was between HB population in 
Lesser Kingan Mountains and HD population in Changbai 
Mountains, while the largest gene flow was between MDJ popu¬ 
lation and JH population (JV e m=12.35) that both of them are be¬ 
longed to the same water system, the Mudanjiang river system 
(Table 2). HD population and JH population are the closest in 
distance, and belonged to Songhua River system (Table 2), 
therefore, gene flow between them was also strong (A e m=12.01). 
Such genetic distance agreed with that inferred from morpho¬ 
logical features such as body length, fore- or hindleg length and 
body mass (Xie et al. 1999). These genetic findings are generally 
consistent with their geographic locations and environmental 
similarity. 

Conservation and management proposition 

Management unit is essential for genetic diversity preservation of 
species. Management unit was usually defined based on the sig¬ 
nificant difference in allele frequency of nuclear DNA and/or 
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mitochondrial DNA, regardless of the occurrence of systemic 
differentiation between populations or between distribution re¬ 
gions (Moritz 1994). The results of this study show significant 
genetic differentiation and weak gene flow between the popula¬ 
tions distributing in Lesser Khingan Mountains and Changbai 
Mountains. Therefore, populations distributed in two regions can 
be proposed to be two independent management units in term of 
protection and genetic management. Moreover, unique alleles 
observed in three populations namely “Tieli”, "Hebei”, and 
“Tonghua” suggested that they have unique genetic characteris¬ 
tics. This is in accordance with morphological and physiological 
traits (Ying et al. 2008). This suggests special attentions should 
be paid to preserve the genetic uniqueness of these three popula¬ 
tions. 
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